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We calculate the number and energy densities of a quasi-2D Bose-Einstein gas constrained within 
a thin region of infinite extent but of finite width S. The BEC critical transition temperature then 
becomes an explicit function of 5. We use this result to construct a model of high-T c superconduc- 
tivity in cuprates with a periodic layered atomic structure. The predicted behavior of the BEC T c 
agrees with recent experimental findings in severely underdoped cuprates. 

PACS numbers: 05.30.Fk, 74.20.-z, 74.72.-h 

INTRODUCTION 

Since the discovery of high-T c superconductivity (HTSC) in cuprates by Bednorz and Muller in 1986 many 
studies to explain the phenomenon have been reported. Its physical origin is not yet clear. Recent measurements 
[H of photoelectron emission spectral intensities from HTSCs suggest that bound electron Cooper pairs (CPs) 
form already at temperatures higher than the critical T c . This finding is consistent with several theoretical efforts 
[H 0] proposing that HTSC originates from a 2D Bose-Einstein condensate (BEC) of CPs pre-existing above 
T c and formed through a BCS-like phonon mechanism. The 2D character of the phase transition is associated 
with the layered structure of cuprates, which in the case of Y ' Ba^Cu-^O-j - y (YBCO) consists of a succession of 
parallel layers perpendicular to the vertical c-axis with a unit cell of height ~ 12A, and the chemical composition 
CuO — BaO — Cu02 — Y — CuC>2 — BaO — CuO. It is generally agreed that the C1LO2 planes, which in the case of YBCO 
are equidistant from the central Y atom with separation ~ 1.5 A, are mainly responsible for the superconductivity 
in cuprates. Contour plots of the charge distribution derived from energy-band-structure calculations for YBCO 
reveal [f| that the SC charge carriers are mainly concentrated within a slab of width 5 ~ 2.15A about the CuOi plane. 

BCS-like theories 0, 0| contemplate a Hamiltonian Bscs containing the kinetic energies of electrons and holes, 
and a pairing interaction arising from phonon-exchange attractions that overwhelm and Coulombic repulsion. As a 
consequence, bound CPs of electrons or holes with antiparallel spins and charge ±2e arise with an energy-momentum 
relation linear at leading order, rather than quadratic, namely £k — £0 + c\hK with c\ = 2vf/^ in 2D with vf the 
• *h , Fermi speed, while K is the CP center-of-mass momentum wavenumber, and £q the familiar weak-coupling energy 
£0 = — 2Hujd exp[— 2/v N ] for K = CPs, where u>d is the Debye frequency, vq the positive BCS electron-phonon 
. coupling constant, and Nq the electron density of states for one spin at the Fermi level. The linear dispersion relation 
is induced Q by the Fermi sea medium so that CPs propagate like free massless composite particles in the Fermi 
sea (whereas in vacuo they would do so quadratically as ft 2 K 2 / Am* if m* is the effective electron mass). Their Bose 
statistical nature allows them to undergo BEC. 



FORMALISM 



We study a model of HTSC in cupratc materials comprising a quasi-2D BEC of excited CPs of energy £k — £q+ciHK 
constrained to propagate within quasi-2D layers of infinite extent in the a\, 0,2 directions, but of finite width 03 = 6 
in the perpendicular direction. We assume that the CP field satisfies periodic boundary conditions (BC) along the a±, 
<Z2 and 03 directions, although the more restrictive Dirichlct or Neumann BCs may be straightforwardly implemented. 
The average number density and the energy per unit volume of the CP field are given by 

n(T) = UdK a( 9{K ] and u{T) = U dK f K) \ K (1) 

where g(K ) is the exact eigenmode distribution of the field defined by g(K) = J2{ n } ^{K—Kn) with = (2ni7r/ai) 2 + 
(2?i27r/a2) 2 + (2n37r/a3) 2 . By introducing Poisson's summation formula it may be shown that g{K) can be expressed 
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in the more tractable way 
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An alternative derivation based on properties of Bessel functions has been formerly employed to study the Casimir 
energy-momentum tensor in rectangular cavities, both at zero and finite temperatures [6j. We then evaluate ([I]) by 
introducing the CP mode distribution J2]) and the CP excitation energy sk = £k — £o — c\hK. The Bose-Einstein 
denominators in the integrals render a rapid convergence so that they may safely be extended to infinity. The integrals 
can be computed by expanding the integrand in powers of ze~ x with the result 
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It is easily checked that the usual thermodynamic limit is attained by considering the terms with mi = m 2 = m 3 = 
in the former expressions. On the other hand, for CPs constrained to move within a layer of finite width 8 <C a\, a 2 
but unconstrained along the infinite ax, a 2 directions we must set m\ = m 2 = in® and @. In that case, 
we introduce the dimcnsionlcss thickness variable r/ = k B TS/hci and the remaining summations over m 3 may be 
performed analytically to give 
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THICK- AND THIN-LAYER LIMITS 



The thick-layer limit given by r) S> 1 represents unconstrained propagation of CPs throughout the entire volume of 
the material as in conventional 3D superconductivity. In this limit h m {rf) — > 1, f m (il) ~ * lj 9m(v) ~ > 1 and we recover 
well-known expressions describing 3D BECs 0, Q : 



n(T) = no(T) + — g- ^ -3 and u(T) = 3— 3- ^ 
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The critical temperature T c follows from the conditions no(T c ) — ► and z(T c ) — > 1 leading to k B T^ D = 
[TT 2 h 3 4n 3D /C(^)] 1/3 where £(n) is Riemann's C-function. The molar heat capacity C(T) = R(nk B )- 1 du(T)/dT 

(with R the gas constant) is straightforwardly obtained from (j6|) as C(T) = ( 12 ( ^ 3 j 4 " > J (t") ^ or ^ < an d 

C(T) = ( 12 m) 4 " > ) (^~) _ ( 12 c(2) 3 ^ ) ^ or T > T c which is consistent with measurements in conventional 3D super- 
conductors [8| since at T = T c the heat capacity shows a discontinuous drop AC = 6.57R indicative of a second-order 
phase transition. 

In the thin-layer limit 77 <C 1 associated with HTSC we get h m (r]) ~ * 0, fm{i]) — mi:/2rj, and g m {r]) — rmv/3Tj. 
Simple algebra leads to 
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where n 2D = nS and u 2D = u8. The critical BEC temperature is now given by ksT 20 = [2nh 2 c\n 2D / X{2)] 1 ^ 2 where 
C(2) = tt 2 /6. In this case the molar heat capacity is C{T) = [6i?C(3)/C(2)] (T/T c ) 2 for T < T c but it must be 
evaluated numerically for T > T c . It turns out that the C(T) is continuous at T = T c although its derivative dC/dT 
is discontinuous. The linear behavior C{T)/T oc T for T < T c is characteristic of cuprate materials 

Of crucial importance in evaluating T c is to reliably estimate the fraction of charge carriers that actually contribute 
to the supercurrcnt. The charge carrier density is usually determined from measurements of London penetration 
depth A a b along the CuC>2 planes. It gives an estimate of the supercurrcnt that causes partial rejection of an applied 
external magnetic field in the superconductor. Within the framework of the present model the supercurrent is due to 
massless-like CPs of charge 2e moving with the CP speed ci, so that the surface supercurrent J s = n 2D (2e)cik with 
k = k/k 0. A straightforward calculation pj shows that n 2D = (e 2 /32irclc 2 )8A 2 ) /Tiu>D^ab where c is the speed of 
light and Ao is the zero-temperature energy gap. The final expression of the BEC critical temperature is 

Tc he / 3,5 N^Aq 
2nk B e \2Tiuj d ) X ab 



CONCLUSIONS 

Introducing in © the YBCO parameters tabulated in Ref.[j| Q D = 410 K , A = 14.5 mcV, and 8 = 2.15 A @ 
we get the relation T c = W.7Q/\ a b (/xm-K) which accurately reproduces experimental data reported by Zuev et al. 
[To| | in measurements performed in YBCO films with T c s ranging from 6 to 507^. They conclude that, within some 
noise their data fall on the same curve X~ 2 oc T 2 - 3±0A regardless of annealing procedure, oxygen content, etc. In an 
independent study, Broun et al. [ll[ found that their samples of high-purity single-crystal YBCO followed also the 
rule T c oc A" 1 . A forthcoming paper Q gives a more detailed discussion of the model presented here applied also to 
several other cuprates including the ^-wave extension of the present formalism valid only for 1 = 0. 
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